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Abstract 

The nature of a jet's fragmentation in heavy-ion collisions has the potential to cast light on 
the mechanism of jet quenching. However the presence of the huge underlying event complicates 
the reconstruction of the jet fragmentation function as a function of the momentum fraction z 
of hadrons in the jet. Here we propose the use of moments of the fragmentation function. These 
quantities appear to be as sensitive to quenching modifications as the fragmentation function 
directly in z. We show that they are amenable to background subtraction using the same jet-area 
based techniques proposed in the past for jet p t 's. Furthermore, complications due to correlations 
between background-fluctuation contributions to the jet's p t and to its particle content are easily 
corrected for. 
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1 Introduction 



In heavy-ion physics one is interested in studying the hot and dense medium that is formed by 
the high energy collisions of heavy nuclei. One way of accessing such properties is to look at 
how the characteristics of a number of 'hard probes' are modified by their passage through this 
medium. While early examples of hard probes consisted of quarkonia or single high transverse- 
momentum (pt) hadrons, the performance and the acceptance of the Large Hadron Collider (LHC) 
detectors (ALICE, ATLAS and CMS) allow for the in-depth study of fully-fledged jets. A number 
of experimental results on jets in heavy ion collisions have already been published following the 
first two heavy ion runs in 2010 and 2011. Of particular interest was the observation by ATLAS [1] 
and CMS [2] |3] of a sizeable asymmetry in the transverse momentum of dijet pairs, which was 
interpreted (modulo some caveats [1] concerning a possible partial spoofing role of background 
fluctuations) as evidence of some amount of quenching, i.e. of a reduction of a jet's momentum due 
to its interaction with the medium. A better understanding, also at the quantitative level, of this 
phenomenon would naturally allow one to better constrain the characteristics of the medium itself. 

After studying an inclusive observable like the quenching of a full jet, the next logical step is 
to observe the momentum distribution of its constituents, i.e. its fragmentation function (FF). One 
would like to ascertain if and how the FF is modified, due to the interactions with the medium, 
with respect to the FF of jets produced in proton-proton collisions and fragmenting instead in the 
vacuum. 

Measurements of jet fragmentation functions have been performed since the first observations 
of hadronic jets. LEP 0E1E] an d Tevatron [8, 9j experiments studied them in detail in e + e~ and 
pp collisions respectively and compared them to accurate calculations in perturbative QCD, finding 
very good agreement. LHC experiments [10} [TT] more recently measured them in pp collisions. 

In a heavy ion context, theoretical discussion as to how the fragmentation function could be 
expected to be modified by the medium were put forward, for instance, in [I~2l [TBI [TH [T5l fl6| [T7 \ [T8] . 
and the question can be examined also with the help of Monte Carlo programs that simulate jet 
fragmentation in a medium [THJ [20], EJJ [22J [23]. The study of jet FFs in a heavy-ion context is 
however complicated by the presence of a large and predominantly soft background due to the 
underlying event produced by the nucleons not directly participating in the hard interaction. The 
background acts in essentially two ways: it adds many soft particles to the jet — an order of 
magnitude more than are normally present in a jet — severely contaminating the fragmentation 
function at low hadron momentum fractions z; and it can bias the measured pt of the jet, even after 
subtracting the expected average pt shift, skewing the value one reconstructs for a given hadron's 
momentum fraction. 

Experimentally, there have so far been measurements of jet FFs by three experiments, in AuAu 
collisions at RHIC by the STAR collaboration |24[ I25j . and more recently in PbPb collisions at 
the LHC by the CMS [11] [26] and ATLAS [271 [28] collaborations. In each case, the measurements 
were limited to hadron pt's above 1 — 4 GeV, where soft contamination is less important. Different 
approaches were used to address the biases in the determination of z, for example unfolding this 
effect in the case of ATLAS, or folding the effect into the pp reference FF in the case of CMS, 
the only experiment which goes down to 1 GeV. The largest changes of heavy-ion FFs relative 
to reference pp FFs appear to be in the region of small z, where uncertainties associated with the 
background subtraction grow. It would clearly be of interest to probe the fragmentation functions at 
lower z and with higher accuracy, especially as other measurements also indicate that soft particles 
appear to play an important role in quenching [2]. 

The purpose of this article is to suggest the use of a new observable that may help provide 
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Figure 1: Jet fragmentation functions (versus z and £) and moments (versus N) in proton-proton 
collisions at the LHC (^/s/vw = 2.76 TeV), obtained without quenching, from Pythia 6, and with 
quenching, from Pyquen. 



an alternative way of understanding fragmentation patterns in heavy-ion jets: momentum-fraction 
weighted moments of the fragmentation functions. 

Moments are often used in theoretical calculations, and have been helpful in the past in di- 
agnosing and resolving experiment-theory discrepancies [29], but to our knowledge have not been 
examined so far in heavy-ion collisions. Among their benefits, one can note that they are help- 
ful in directly relating jet spectra with hadron spectra and also that they are sensible objects on 
a jet-by-jet basis (whereas event-by-event, fragmentation functions are simply sums of individual 
delta- functions, one for each hadron in the jet). The latter point makes it feasible to measure 
correlations within an event between any given moment of the background fragmentation spectrum 
and fluctuations in the background's pt- As we shall see, it seems that this can be of significant 
benefit in extracting the soft part of the jet's fragmentation and also in reducing systematics in the 
extraction of the hard part. 



2 Representations of the fragmentation function 

A jet fragmentation function can be defined as the distribution dN^/dz of the momentum fraction 

Pt,h n , 

z = , (1 

Pt 

of hadrons in the jet, where pt t h is the transverse momentum of the hadron and pt is the transverse 
momentum of the jet0 

Figure Q] shows, in its left-most plot, the typical shape of the z distribution (normalised in this 
case to the total number of jets Aj e t used in the analysis), given for anti-kt (R = 0.4) jets [30] 
with pt > 100 GeV in pp collisions at a centre-of-mass energy of 2.76 TeV. The middle plot shows 
identical information, represented as a differential distribution in £ = \nl/z. The z representation 
helps visualise the hard region of the FF, while £ devotes more space to the soft part. Two curves 
are shown: one of them, solid blue, labelled "Pythia 6", is a pp reference curve obtained from 
Pythia 6.425 |31] with its virtuality-ordered shower. The other, dashed green, was obtained with 
the Pyquen program [19, 32j, which modifies Pythia showering so as to simulate quenching. We 



1 Various other definitions are also in use, for example replacing transverse momenta with 3-momenta and taking 
the projection of the hadron momenta along the axis of the jet. 
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have used it with settings corresponding to — 10% centrality. Its effect on the FF is of the same 
order of magnitude as the effects seen experimentally |26t [28] and it therefore provides a useful 
reference when establishing whether FFs are being reconstructed with sufficient accuracy by some 
given procedure. 

The right-most plot of Fig. Q] shows moments of the fragmentation functions. The -/V th moment, 
Mjv, of the fragmentation function is given by the integral 

\. N dN h _ 1 [°°_ N6 dN h 



iVjet J dz iV jet Jo V 7 



In practice, the moments for a single jet can be calculated as 



M jet = ^iejetPt,, ^ (3) 



where the sum runs over all the jet's constituents. The results can then be averaged over many jets, 
so that Mjy = (-^v^jets- Obviously, Mo represents the average particle multiplicity in a jet, and M\ 
is equal to one by virtue of momentum conservation if one measures all hadrons, as we assume here 
(taking 7ro's to be stable) H If instead only charged tracks are used in the numerator, then it is clear 
that Mi will be significantly below 1. There is another value of N that is of special interest: given 
a jet spectrum da^t/dpt that falls as Pt~ n , the ratio of the inclusive hadron spectrum and inclusive 
jet spectrum is given by M n -\. Thus, M^ l /M pp _ l corresponds to the ratio of (charged-) hadron 
and jet Raa values (in the approximation that n is exactly independent of pt), 



M n-1 _ R AA 
±v± n— 1 -"-AA 



(5) 



For p t in the range 100 — 200 GeV, at \/s NN = 2.76 TeV, the relevant n value has some dependence 
on pt and is in the range n = 6 — 7, corresponding to N = 5 — 6. 

In representing the moments in figure QJright), we include a factor ((N + l)/2) 4 , which allows 
a broad range of N values to be shown on a linear vertical scale. The same features visible in the 
plots versus z and £ are visible versus N too, for example that Pyquen leads to higher multiplicities 
than Pythia at small N (corresponding to z < 0.05 or £ > 3) and large N (z > 0.5, £ < 0.7), and a 
slightly reduced multiplicity at intermediate N (z ~ 0.2, £ ~ 1.5). 

To help understand the quantitative relationship between N and £, one may examine figure El 
a colour-map that shows as a function of ./V the contribution to the Mn moment from each £ value. 
It shows clearly how large £ values dominate for low (and vice- versa). This £,N relationship 
depends to some extent on the shape of the fragmentation function and it is given for the same 
Pythia 6 fragmentation function that was used in Fig. [TJ 



2 The relation Mi = 1 is exact only if one defines the transverse momentum of a jet as the scalar sum of the 
transverse momenta of the constituents, so that the moments are given by 

v-* JV 

L.iPt,i (4) 

This is the choice that has been made in this paper. Using the transverse component of the jet momentum rather 
than the scalar sum leads to small violations (of the order of a fraction of one percent, for the jet radius R — 0.4 
used in this work) of the Mi = 1 relation. If the fragmentation moments are measured with only charged hadrons, 
we will still assume that the denominator of Eq. Q is determined with all particles in the jet or, equivalently, with 
all calorimeter towers. 
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Figure 2: Representation of the £ values that contribute dominantly to the Mjv integral for a given 
N, shown as a function of N. Shown for the Pythia 6 results and cuts of Fig. Q3 



3 Simulation and reconstruction tools 



Before proceeding with the heavy-ion analysis, it is useful to detail the simulation tools that we 
use. 

What we call 'hard' jets in QCD are simulated in proton-proton (pp) collisions using Pythia 6.425 [31] 
in dijet mode with the DW tune [33jH We have not included any pp underlying event (UE), its effect 
being minimal in this context anyway. Jets including quenching effects have been generated using 
Pyquen [HI [32], v 1.5.1, as included in Hydjet [HUES]. We consider all hadrons, not just charged 
tracks, and take 7Tn's to be stable, so that we are still considering genuine hadron distributions. 

The heavy-ion background (which we will also call heavy-ion underlying event) is simulated 
using Hydjet vl.6 for 0-10% centrality and, where needed, it is superimposed to the hard event 
generated by PythiaEI The jets observed in this combined event will be denoted as 'full' jets. 

All events are generated for either pp or lead- lead (PbPb) collisions at the LHC, with a centre- 
of-mass energy of 2.76 TeV per nucleon-nucleon collision. 

Jets are reconstructed using the anti-fcj algorithm [30] with R = 0.4, as implemented in Fast- 
Jet [HU [37]. In estimating the jet's pt, the HI background is subtracted using the median/area 
based techniques introduced in [38l [39j 130], which we implement with kt jets [H] (R = 0.4) to 
estimate p, the background transverse momentum per unit area. At most the two hardest jets 
passing a hard cut on the subtracted transverse momentum (see below) are subsequently used for 
the fragmentation function analysis. The cuts that we shall use are 100 GeV and 200 GeV. 



4 Impact of HI background and its subtraction 

As explained in the Introduction, the addition of the heavy-ion background has the potential to 
modify a measured fragmentation function in two ways. Firstly, the jet's pt is modified, affecting 
the normalisation of z in Eq. ([1]). Secondly, the heavy-ion background adds many extra particles 
to the jet, predominantly at low momenta. 

3 This tune is not the most up to date; however it is not unrealistic for the LHC and has the characteristic that it 
is based on Pythia's virtuality ordered shower, which is a prerequisite for use with Pyquen. 

4 In their notes, CMS and ATLAS often call this procedure 'embedding' or 'overlying' respectively. 
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To appreciate the impact of the extra particles in the jet from the heavy-ion background, it is 
instructive to first examine the "Pythia+Hydjet" dashed red curves of Fig. [3l These show the FF 
extracted in heavy- ion collisions, without any subtraction of the background contribution to the 
FF, but always using a z value defined such that the jet's pt has been corrected for the expected 
HI background contamination (as is standard in the experimental measurements). At this stage we 
will not perform any unfolding to account for fluctuations in the HI background. 

One sees how the FF acquires a 'bump' in the soft region, which lies at larger £ (smaller z) than 
the maximum in the original pp result (blue solid line) and is up to two orders of magnitude higher. 
For the moments, this presence of the background is seen as a steep increase in the small- N region, 
taking the curve far off the scale. 

At large z (small £), the impact of the addition of the background is only barely visible, as 
might be expected given that it is dominated by soft particles. However, this visually small effect is 
partially an artefact of the logarithmic scale used to show the FF versus z. Considering instead the 
moments, one sees that there is a non-negligible reduction in the FF at large N. This is perhaps 
surprising given that the background adds particles. It is a feature related to the interplay between 
background fluctuations and the steeply falling jet spectrum. It is well known by the experiments 
and we will return to discuss it in section [5j 

4.1 z-space subtraction 

The traditional approach to background subtraction from a jet FF involves the construction of a 
distribution (in z or £ space) intended to approximate that of background-only particles, and the 
subtraction of this distribution from the measured one. The way the background-only distribution 
is determined can vary, the simplest one probably being to consider a region of the event that is 
expected to be little affected by the hard jets, and measure it there. 

To illustrate z-space subtraction here, we measure the distribution of hadrons in two regions 
transverse in azimuth with respect to the axis defined by the dijet event. Event by event, and jet 
by jet, we subtract those distributions (measured in a patch of phase space with an area equal to 
that of the jet that one is considering) from the jet fragmentation function. While exact experi- 
mental procedures differ in the details, most choices lead to similar results here. Perhaps the main 
distinction of the experimental procedures is that they sometimes address issues related to flow, 
which for simplicity we neglect here in our z-space subtraction. 

The results of this subtraction are shown in figure[3]as green open circles, labelled 'z-subtracted' 
(the moments of these ^-subtracted results are also shown) . One sees how these curves come closer to 
the pp results at small z and small N than do the unsubtr acted (dashed red) ones. The agreement 
improves as the jet pt threshold is increased. However, even with jets of pt > 200 GeV, this 
procedure still falls short of an accurate reconstruction of the hard FF. In fact, the region in z 
where the subtracted FF works well barely extends beyond that selected by simply truncating the 
unsubtracted FF at low pt t h so as to avoid the region dominated by the background. 

4.2 iV-space subtraction 

An alternative approach is to directly subtract the moments Mjy of the fragmentation function. 
This can be done by extending the jet-area based techniques for background estimation introduced 
in [38] and implemented in [36]. 

Let us first recall the procedure for subtracting the background contribution to a jet's pt- One 
first determines the background pt density per unit area, p, for the event (or just part of the event 
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Figure 3: Jet fragmentation functions shown for plain Pythia, with the addition of the heavy-ion 
background (Pythia+Hydjet) and after subtraction of the heavy-ion background (z-subtracted and 
iV-subtracted). For the results including the heavy- ion background, the jet p t used to define z is 
always that after subtraction of the heavy-ion background. As in Fig. [IJ we show the results as a 
function of z (left), £ (middle) and for the moments versus N (right). The upper (lower) row has 
a jet pt threshold of 100 GeV (200 GeV). 



in the vicinity of the jet). One accomplishes this by dividing the event into patches of similar 
rapidity-azimuth area, e.g. by running the inclusive kt algorithm and taking all the jets it finds. 
Then p is obtained by taking the median across all patches of the ratio of pt to area (A patc h) for 
each patch: 

j. I Pi, patch I fa \ 

p = median < — p- > . (t>J 

patches [ ^4patch J 

The median serves to limit biases from any hard jets among the patches. A given jet's full transverse 
momentum, pt -Mi, which includes genuine jet particles and background particles, is then corrected 
for the background particles by subtracting an amount given by the product of p and the jet's area 
A: 

Pttn = Pt,™ ~ P A • (7) 
We use patches in an annulus (or 'doughnut') of outer radius 3R and inner radius R, centred on 
the jet of interest. This choice helps to limit corrections due to flow and accounts for the rapidity 
dependence of the background!! 

For the purpose of correcting FF moments, the procedure we propose is quite similar. One first 
determines the expected background contribution per unit area to a given moment (or rather to its 

Alternatively one could determine a global p for the event and modulate it for rapidity and azimuthal dependence, 
as discussed in (35], and similarly to what is done for example by ATLAS |42] , 
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numerator in Eq. ©): 

J ^-^iGpatch Pt.i I z„\ 

Pn = median < j- > . (8) 

patches ^4patch J 

For example, for N = this will be the median particle multiplicity per unit area. The subtracted 
FF moment is then obtained by separately taking the numerator and denominator of Eq. ([3]), as 
measured in the full event, and respectively subtracting p^ A and pA: 

A/f sub _ EigjctPM ~PnA 

(Pt.fuii - pA) N 

The results for M^ ub are shown in figure [3] in the right-most plots, orange crosses labelled l N- 
subtracted'. One can see that they are neither better nor worse than the corresponding z-subtracted 
ones, which have been translated to iV-space and drawn in the same plots for direct comparison 
(green circles). 

In order to provide a fragmentation function reconstruction that is better than that given by 
the standard z-space method, we introduce in the next section an improved background subtraction 
method that is most straightforwardly applied in moment space. 



5 Improved background-subtracted fragmentation function 

The discrepancies that we have observed can, we believe, largely be attributed to fluctuations of 
the background. 

Background fluctuations mean that neither the jet's pt nor the numerator of the FF moment are 
perfectly reconstructed for any given jet. When selecting jets above some pt threshold for a process 
with a steeply falling jet spectrum, it is favourable to select jets slightly below the pt threshold, 
but which have an upwards background fluctuation (cf. the discussion of section 5 of [10]). One 
consequence of this is that the denominator in Eq. ([9]) tends to be larger than the actual jet pt- The 
larger the value of N, the greater the impact of this effect. One can also understand it as resulting 
in an underestimate of the z fraction, leading to the FF being spuriously shifted to lower z. Such an 
effect is already well known: CMS for example, when comparing with pp FFs, explicitly applies 
a correction to the pp FFs to account for the smearing of the denominator of the z variable that is 
expected when the jet's momentum is reconstructed in a HI environment. ATLAS [28J carries out 
an unfolding to account for this. 

A second consequence of fluctuations affects mostly the low-z, or low-A^ region of the FF: the 
upwards fluctuations of the background pt that cause a jet to pass the pt cut even when it is below 
threshold also tend to be associated with upwards fluctuations of the multiplicity of soft background 
particles. It is this effect that causes the FF to be overestimated for low values of ivjj 

One way of verifying the above interpretation is to consider 7+jet events, selecting events based 
on the photon pt and normalising the fragmentation function z also to the photon pt. In this case, 
the background-induced fluctuations of the jet's reconstructed pt are of no relevance, and in explicit 
simulations, we have found that plain subtraction is already quite effective. 

For the dijet case, in the limit where the background fluctuations are reasonably small compared 
to the transverse momentum of the jet, it is possible to devise a simple correction in moment space 

6 As this letter was being finalised, ATLAS presented preliminary results [28] in which they correct for this by 
multiplying the expected background contribution by a ^-independent factor that is a function of the jet pt- 
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for both kinds of fluctuation-induced bias. In order to do so, we start by rewriting the subtracted 
moments of the jet FF in Eq. ([9]) as 



sub _ Y.iPt,i-pNA _ s N 
MN - (Pt M i-pA) N ~ S? (10) 

where we have introduced the shorthands Sn for the subtracted numerator and Si for p^n, as 
defined in Eq. (|7|). Given our assumption that background fluctuations are moderate, we can locally 
approximate the hard jet cross-section by an exponential rather than by a power in order to facilitate 
our analytical working, 

H( Pt ) = ^ = ^eM-Pt/^)- (11) 
apt p 

Next we take a Gaussian approximation for the spectrum of background transverse-momentum 
fluctuations from one jet to the next. Denoting the fluctuation by qt, the distribution of qt is then 

, . dP 1 f of \ 

where a is a parameter that describes the size of fluctuations from one patch of area 1 to another. 
It can be extracted directly from the event in a manner similar to p [38, 43J. Its value is of the order 
of a ~ 18 GeV for our PbPb LHC (^Snn = 2.76 TeV) simulations (the same simulation leads to 
fluctuations of <Tj e t — 11 GeV for anti-fct R = 0.4 charged-track jets, in good agreement with the 
measurement from ALICE |44j ) . 

We further introduce the variable Qn to denote the difference between the actual background 
contribution to Sn in a specific jet and the expected contribution, pnA, i.e. 



Qn=\ V k»\- PN A, (13) 




where the sum runs just over the background constituents fc^ of the jet. By construction, Qi = qt- 
In practice Qn cannot be determined for a single jet, since we don't know which particles are 
the background ones, but its statistical properties can be determined by looking at many jets. 
The fluctuations of Qjv are not independent of the fluctuations q t of the background's transverse 
momentum: there is a correlation coefficient rjy between them, defined as 

Cov (g t , Qn) n ., 
r N = ; , (14) 

VVar(^)Var(Q^) 

where Var(X) is the variance of the variable X and Cov(X,Y) the covariance of X and Y. Using 
the fact that Var((ft) = a 2 A, Var(QAr) = cr N A, we have that the average value for Qn as a function 
of q t is 

/n \( \ C(w(q t ,Q N ) o-JV / 1K s 

(QN){qt) = — 77 , s — qt = rjv — qt ■ (15) 
Var(g t ) a 

With these ingredients we can now correct for the fluctuation effects as follows. Firstly, we note 
that the cross section for a given reconstructed Si (i.e. subtracted jet pt) is 

J| = J dq t H(Si - qt)B(qt) , (16) 
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where the integral allows us to deduce the probability distribution for the actual fluctuations q t 
given S\: 

dP 



dq, 



H(S 1 - q t )B(q t ) 



Si fd q >H(Si-q' t )B( q > t ) 



corresponding, with our approximations for H(pt) and B(qt) in Eqs. (jlip and f)12|) respectively, to 
an average qt of 

rr 2 A 

(it) = — • (18) 
A* 

If we measure a certain value Sn in a jet, then as a function of the qt fluctuation, the expected true 
hard contribution to it is 

S%* rd = S N - (Q N )(qt) = S N - r N ^q t , (19) 

a 

where we have averaged over possible Qn values, given the qt fluctuation. To obtain our estimate 
for M]jr ard , this should be normalised by the ./V th power of the true hard pt of the jet, (Si — qt) N - 

One subtlety here is that this is an estimate for M]^ ard in hard jets with true pt = S\ — qt- However 
because M/v is a slowly varying function of pt, by taking the result as contributing to M^ rd (Si) 
rather than M^ rd (Si — qt) we make only a small mistake, of the same order as other terms that 
we shall neglect 

Retaining the terms linear in qt in Eq. (I20j) and averaging now over possible qt values, making 
use of Eq. (fT8|). leads to the following prescription for an "improved" subtracted Mn(Si), corrected 
for fluctuations effects up to first order in qt/S\: 



M su bMP = M ^x[l + N ^ ) _ r N ^f . (21) 



,a 2 A\ aa N A 
Shi J /xSf 

This is simpler than the corresponding correction would be directly in z space, in particular because 
in z space the correction to one bin of the fragmentation function depends in a non-trivial way on 
the contents of nearby bins. One can think of the advantages of moment space as being that the 
correction to a given N value does not depend on M^r at all other values of N, and that it is 
straightforward to account for correlations between fluctuations in the jet-pt and in the moments. 

Note that in a real experimental context, calorimeter fluctuations of the reconstructed jet and 
background pts would have an effect akin to increasing a and decreasing the correlation coefficient 
rfq. The in-situ methods that we use for the determination of a, <jn and would automatically 
take this into account. Noise-reduction methods in the reconstruction of the jet pt, as used by 
CMS [2], would have the effect of reducing a (and probably also r^). However noise reduction is 



7 This can be seen by observing that Mn satisfies a DGLAP-style equation for its evolution dAfjv / 'dlnp t ~ a s Mjv. 
Given that qt is itself small compared to pt, Mn for jets with p t — Si — qt differs from that for jets with p t = Si 
by a relative amount ~ a s qt/Si, which we can neglect in the same way that we neglect O (qt) terms. If we wanted 
to improve on this approximation, then one approach might be, for each event, to use the middle expression in 
Eq. Q20p and assign it to Mn at p t = Si — (qt)- We have not, however, investigated this option in detail and other 
improvements would probably also be necessary at a similar accuracy, e.g. taking into account deviations from the 
simple exponential and Gaussian approximations that we have used for H(p t ) and B(qt). 
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Figure 4: The quantities pn, <7jv and the correlation coefficient r^r, shown as a function of N for 
— 10% central Pb Pb collisions -^/snn = 2.76 TeV as obtained from simulations with Hydjet. 



likely to complicate the meaningful determination of r^r, since it acts differently on pure background 
jets as compared to jets with a hard fragmenting component. 

The correction in Eq. ([2T|) can be applied jet-by-jet to correct for the fluctuation effects. It 
requires the prior knowledge of the slope \x of the jet cross-section, which can be obtained from pp 
data, or from simulationsO All the other ingredients that enter this equation (a, A, Sn, <tjv, rjv) 
can instead be determined event-by event or jet-by-jet. In practice we determine a, and from 
the ensemble of jets contained in an annulus (or "doughnut" ) of outer radius 3R and inner radius 
R, centred on the jet of interest. Typical values of pjy, <tjv and rjy are presented as a function of N 
in Fig. @J 

We show in figure [5] the result of applying Eq. (|21jl to our subtraction in moment space, for two 
jet pt thresholds. The solid blue curve {pp reference) and orange crosses (N-subtracted) are identical 
to the results in the rightmost plots in figure El In addition, figure [5] also displays results obtained 
using Eq. (|2ip . shown as red diamonds. One sees how the quality of the agreement with the 'hard' 
blue curve is markedly improved. At low N it is the last, additive, term in Eq. (|2ip that dominates 
this improvement, accounting for the correlation between a jet's reconstructed pt fluctuations and 
the fluctuations in the background's contribution to the moment; at high N it is the multiplicative 
Na 2 A/S\p term that dominates, correcting for the fact that the more common upwards fluctuations 
in the jet's reconstructed pt cause the fragmentation z value to be underestimated. 

Figure [5] also shows (green dashed curve) the jet fragmentation function predicted by the quench- 
ing model used in Pyquen. One observes that the remaining deficiencies of the reconstruction are 
significantly smaller than the difference between the unquenched (solid blue) and the quenched 
(dashed green) FFs, pointing to a potential discriminating power. This is to be contrasted with 
subtraction without improvement (orange crosses), which especially in the soft region, N < 1, fails 
to describe the blue curve sufficiently well to tell whether quenching (as predicted by Pyquen) is 
present or simply that an imperfect reconstruction is taking place. This serves as an illustration 
that the improved subtraction may now be sufficiently good to allow one to discriminate a quenched 
FF from an unquenched one. 

We conclude this section by noting that an implementation of the tools needed to implement the 

8 In practice, fi depends on p t and should be taken at the scale Si — qt in the integrand. However, /i varies slowly 
with p t and can easily be taken at the fixed scale pt in our small-fluctuations limit. In our analysis /i ranged from 
~ 9 GeV at p t ~ 50 GeV to ~ 28 GeV at p t ~ 200 GeV. At p t ~ 100 GeV we had n ~ 16 GeV. 
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Figure 5: Jet fragmentation function moments, showing the plain Pythia result, the result after 
embedding in Hydjet and applying plain subtraction moment-space subtraction ( "N-subtracted" ) 
and after the additional improvement to account for correlations ( "+ correl" ) , Eq. (|2ip . A quenched 
result ("Pyquen") is also shown, to help give an indication of the order of magnitude of quenching 
effects as compared to residual misreconstruction effects. 



jet fragmentation function subtraction in moments space, as well as the fluctuations unfolding im- 
provement, will be made available as a Fast Jet add-on from |http: //fast jet . hepforge . org/contrib, 

6 Conclusions 

In this paper we have suggested that it could be of interest to study jet fragmentation-function 
moments in heavy-ion collisions. Individual moments contain specific information such as jet- 
hadron multiplicities, the momentum fraction carried by charged particle and can also provide 
insight into the relation between inclusive hadron and inclusive spectra. 

As we have seen, it is quite straightforward to correct moments for effects of fluctuating back- 
grounds, arguably more so than for fragmentation functions expressed in terms of momentum 
fractions. The correction procedure we discussed is amenable to systematic improvement, for ex- 
ample by accounting more completely for the shape of the jet pt spectrum. Such extensions are 
best considered as part of a more detailed analysis, for example in the context of a full experimental 
study. 
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